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ABSTRACT 

We present a universal linear correlation between the stellar mass and surface brightness (SB) 
of galaxies at 0.3 < z < 3, using a deep if-band selected catalog in the GOODS-North region. 
The correlation has a nearly constant slope, independent of redshift and color of galaxies in 
the rest-z frame. Considering unresolved compact galaxies, the tight correlation gives a lower 
boundary of SB for a given stellar mass; lower SB galaxies are prohibited over the boundary. 
The universal slope suggests that the stellar mass in galaxies was build up over their cosmic 
histories in a similar manner irrelevant to galaxy mass, as oppose to the scenario that massive 
galaxies mainly accumulated their stellar mass by major merging. In contrast, SB shows a strong 
dependence on redshift for a given stellar mass. It evolves as ~ (1 + 2;)~2 0~-o.8^ jj-^ addition to 
dimming as (1 + z)^ by the cosmological expansion effect. The brightening depends on galaxy 
color and stellar mass. The blue population (rest-frame U — V < 0), which is dominated by 
young and star- forming galaxies, evolves as ^ (1 + z)^'^-^^'^-^ in the rest-V^ band. On the other 
hand, the red population {U — V > 0) and the massive galaxies (M, > lO^^M©) shows stronger 
brightening, z)~^'^^^''^ . Based on the comparison with galaxy evolution models, we find that 
the phenomena are well explained by the pure luminosity evolution of galaxies out to z ^ 3. 

Subject headings: galaxies: evolutions — galaxies: fundamental parameters (surface brightness) — galax- 
ies: high-redshift — infrared: galaxies 

1. Introduction forming and merging histories of galaxies. The 

surface brightness (SB) of galaxies is also key to 

Galaxy evolution from high-z to the present studying galaxy structure. For example, Kor- 

time is one of the central issues to be understood ^-^^^^^ (^977) showed a tight correlation between 

in ACDM paradigm. Among many parameters to gg ^^dius parameters for compact and normal 

characterize galaxies, stellar mass, which is accu- galaxies. Since then, many authors have studied 

mulated in a galaxy with aging, is one of the most ^^e evolution of SB in the high-z universe (e.g., 

fundamental quantities for understanding the star- Jorgensen et al. 1995; Roche et al. 1998; Lilly 
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et al. 1998; La Barbara ct al. 2004; Barden 
et al. 2005; Ferreras et al. 2009; Trujillo 2004, 
2006; Cimatti et al. 2008; van Dokkum et al. 2008; 
Muzzin ct al. 2009: Damjanov ct al. 2009). The 
SB and its correlation with other structural pa- 
rameters of local galaxies were investigated with 
SDSS galaxies (e.g., Shcn ct al. 2003: Bcrnardi ct 
al. 2003). A direct evidence of SB dimming as 
(1 + z)'^ also intrigues us as a cosmological expan- 
sion test (Tolman 1930; Lubin & Sandage 2001). 

By analogy with the Kormendy relation, the 
correlation of stellar mass with SB of galaxies 
(hereafter called the stellar-mass Kormendy rela- 
tion) in the high-2: universe will give us a clue to 
understanding the history of stellar mass built up 
in galaxies over cosmic time. In that context, we 
will study the correlation and its evolution from 
^ 3 to 2 ~ 0.3, using a deep ii'-selected galaxy 
catalog. 

In §2, we describe the catalog we used. The 
data analysis and the result for the stellar-mass 
Kormendy relation and its evolution are detailed 
in §3 and §4. The results are discussed in §5. We 
use the AB magnitude system, except for Vega 
U — V color. 

2. Data 

We use the if -band selected catalog of the 
MOIRCS Deep Survey (MODS) in the GOODS- 
North region (Kajisawa ct al. 2009, hereafter 
K09), which is based on our imaging observa- 
tions in JHKs bands with MOIRCS (Suzuki et 
al. 2008) and archived data. Four MOIRCS 
pointings cover 70% of the GOODS-North region 
(103 arcmin^, hereafter referred as 'wide' field). 
One of the four pointings, which includes HDF-N 
(Williams et al. 1996), is the ultra-deep field of 
MODS (28 arcmin^, 'deep' field). The deep imag- 
ing gives us SB data in a wide dynamic range, so 
that it is also advantageous for studying low SB 
galaxies at high redshift, where SB dimming due 
to the cosmological expansion is significant. 

Although the detection completeness (90%) of 
the catalog is ii" ~ 25 mag for the wide field and 
~ 26 mag for the deep field, we include fainter 
galaxies, because the if-band limit does not con- 
cern the present study. In fact, the sample with 
K = 25 limit gives the same result within error. 
Instead, we will take account of the size and SB 



limits of galaxies. The FWHMs of the final stacked 
images are 0".46 for the deep image and 0".53- 
0".60 for the wide image. The numbers of galaxies 
are 3717 and 6265, respectively. 

To obtain stellar-mass functions of MODS sam- 
ples, K09 performed SED fitting of the multi- 
band photometry {UBVizJHK, 3.6 /xm, 4.5 fim, 
and 5.8 /xm) with population synthesis models. 
We adopt the results with GALAXEV templates 
(Bruzual & Chariot 2003) and the Salpcter initial 
mass function (see K09 for more details). In the 
present analysis, the near-infrared data (3.6/xm) of 
Spitzer/IRAC are available for most of the sam- 
ple galaxies (96%), so that SED fitting is reason- 
ably reliable for the photometry at rest-z (Aes = 
0.9/im) or shorter wavelengths out to redshift ~ 3. 

3. Analyses 

3.1. Stellar-mass Kormendy Relation 

The surface brightness hk (mag arcsec"^) in 
K band of the sample galaxies is defined as the 
brightness in two times the Kron radius, i.e., 

liK = mK + '2.Mog{TTR^), (1) 

where tuk is total magnitude and R is two times 
the Kron radius r (R — 2r) in units of arcsec. For 
niK, we use MAG_AUTO obtained by SExtractor 
(Bertin & Arnouts 1996), which gives the most 
precise estimate of the total magnitude of galaxies. 
Kron radius is defined as the first moment. 



where r' is isophotal radius and / is fiux at r'. 
About 90% or more of the total light is within R, 
irrespective of galaxy shapes or redshifts (Kron 
1980). 

The half-light radius and central SB of galaxies 
have been investigated in many previous studies. 
However, in contrast, it is less reliable at higher 
redshifts in the present sample because of seeing 
effects. Therefore, we study the total SB in R 
instead. Figure 1 shows fj,K versus stellar mass 
(M^) of all galaxies broken into redshift bins. 

3.2. Point-like Sources 

In order to examine the error of SB, we first 
study the reliability of the Kron radius of galax- 
ies. Figure 2 shows the Kron radius derived by 
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SExtractor versus rriK for all MODS samples. We 
draw a boundary line above the stars, which are 
spectroscopically confirnied. Since spectroscopic 
data were not available for stars fainter than ^ 24 
mag in the present region, we examined the Kron 
radius using artificial stars. The point sources 
were convolved with a gaussian point spread func- 
tion (FWHM=0".5), and then buried in the image 
with the same noise (cr = 30.9 mag pixel"^) with 
that of the deep image. SExtractor was applied to 
the artifacts in the same manner as for the MODS 
catalog. The results are plotted over Figure 2. 
The objects below the boundary in Figure 2 are 
considered as point-like or unresolved sources in 
the following. 

3.3. Monte Carlo Study v^fith Mock Galax- 
ies 

Seeing and background noise strongly affect the 
observed Kron radius, especially for faint or small 
galaxies. To examine the effect, we generated 
mock galaxies with a 1/ 4- law or exponential pro- 
file and buried them in the simulated noise image. 
The galaxies were randomly generated with vari- 
ous magnitudes and effective or scale lengths, then 
analyzed with SExtractor in the same manner as 
for the MODS catalog. The result is shown in Fig- 
ure 3. The observed Kron radii r are seriously af- 
fected by seeing for galaxies with an original Kron 
radius ro < 0".4. The observed SB of the mock 
galaxies with ro < 0".4 is systematically fainter 
than the originals. We note that faint SB galaxies 
are strongly affected by background noise. For ex- 
ample, the large scatter at log{r/ro) < in Figure 
3 is mainly due to low SB galaxies with SB> 26 
mag arcsec"^. Figure 3 also shows that the SB 
of galaxies with observed SB> 26 mag arcsec^'^ is 
unreliable. Therefore, we will exclude the galaxies 
with r < 0".4 and/or SB> 26 mag arcsec"^ from 
the following discussion. 

3.4. K — Corrections 

Although a universal correlation independent 
on redshift is suggested in Figure 1, when we dis- 
cuss the evolution of the correlation as a function 
of redshift, fc-correction should be applied to the 
observed rriK to compare them at the same rest- 
frame wavelengths. The SB, /ii.est(A) in a rest 
frame at wavelength A, is obtained from the ob- 
served data and the model SED, using the equa- 
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Fig. 1. — Correlation of surface brightness in K 
band {^ik) with stellar mass (M,) of galaxies for 
each redshift bin. The galaxies are classified ac- 
cording to the color (rest-frame U — V) and the 
regions (wide and deep) as shown in the bottom 
right frame. 
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Fig. 2. — Kron radius r and total magnitude rriK 
for the MODS catalog. Squares show spectroscop- 
ically confirmed stars. The results of SExtractor 
for artificial stars are plotted by filled circles. The 
solid line delineates the upper boundary of point- 
like or unresolved sources. 
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Fig. 3. — (a) Kron radius and (b) surface bright- 
ness for mock galaxies. The abscissa is the origi- 
nal values and the ordinate is those derived with 
SExtractor. Crosses (blue) and circles (red) show 
galaxies with 1/4-law and exponential profiles, re- 
spectively. The galaxies with fi > 26 are depicted 
with filled circles and pluses. The solid curve in (a) 
is the observed Kron radius r — 0.4. The galaxies 
below the curve are noted with filled triangles and 
asterisks in (a) and (b). The galaxies below the 
solid curve at r = 0.4 in (a) and above the dashed 
line at /i = 26 are excluded from the present study. 
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Fig. 5. — Same as Figure 4, but for the rest-z 
frame. 
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Fig. 4. — Same as Figure 1, but for the rest-T^ 
frame obtained with Equation 3. Typical errors 
are depicted by crosses for mass bins. Solid, dash, 
and dash-dot lines are the linear fits for all, blue 
([/ - y < 0), and red {U - V > 0) samples, re- 
spectively. 
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Fig. 6. — Same as Figure 4, but for discarded 
galaxies. The same regression lines as in Figure 
4 are shown for reference. 
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tion, 



with a linear regression, 



fXreA^)=mK + k{X) + 2.5log{TTR^)-2.5log{l + zf, (3) /x.est = a^offM* - 10) + /x.^^t, (4) 



where fc(A) is the fc-correction defined as the color 
k{\) = (rcst-y — K) or (rcst-z — K) of the best 
SED model fitted to the observations by the 
analysis (K09). 

We apply the fc-corrcction to m.K for each sam- 
ple galaxy to obtain the SB in the rest-^ and z 
bands. The last term (1 + z)^ is the correction 
for the dimming effect due to the cosmological ex- 
pansion, after m.K in AB system is reinstated by a 
factor of (l+z). Figures 4 and 5 shows the correla- 
tion of Mh, and SB in the rest-F and z frames, re- 
spectively. The point-like galaxies (r < 0".4) and 
the galaxies with hk > 26 mag arcsec"^ are dis- 
carded from the analysis. The discarded objects 
are depicted separately in Figure 6 for compari- 
son. Since the sizes of galaxies are suggestive of 
independence on wavelengths in the local universe 
(e.g., Shen et al. 2003; Bernardi et al. 2003) and 
at high redshift (e.g., Cassata et al 2009 at ^ ~ 2), 
no fc-correction for galaxy size is applied in the 
present analysis. 

3.5. Error Estimate 

The errors in stellar mass and fc-corrcction 
mainly originated from SED fitting and the photo- 
metric errors of the observations. In the course of 

fitting to obtain a best SED model with various 
parameters (e.g., star- formation time scale, photo- 
metric redshift for galaxies with no spectroscopic 
redshift available, age, extinction, and metallic- 
ity), the probability distributions of stellar mass 
and fc-correction can be calculated, where the pho- 
tometry and photometric-redshift errors arc in- 
cluded. The error for the observed SB is exam- 
ined with mock galaxies in Figure 3 as a function 
of Kron radius. The photo-z error for galaxies 
with no spectroscopic redshift was studied in K09. 
Prom these errors, we obtained 1-a errors for stel- 
lar mass and surface brightness for each galaxy in 
the rest-y and rest-z bands. The average errors 
for the stellar-mass bins are depicted in Figures 4 
and 5. 

4. Results 

We obtain the least square fit of rest-F (rest-^) 
SB to M* for the galaxies in Figure 4 (Figure 5) 



where nl^J^^. is SB at M* = lO^^M© in the rest 
frame. 

Since our image quality is not high enough for 
classifying the galaxies into morphological classes 
at high-z, we divided the samples into two groups 
according to rest-frame U — V color and obtained 
again the linear fit for each group. Rest-frame 
U — V = corresponds to the color of AO stars, 
so that the bluer group is supposed to be dom- 
inated by young and star-forming galaxies. The 
results of the regression analysis with mean error 
are summarized in Table 1. The dispersion, a, 
of the linear fit is listed in the last column. (For 
galaxies with M* > 10^^ Mq, the average SBs are 
obtained instead of linear fitting.) Figure 7 shows 
a, the slope of Equation 4, as a function of red- 
shift for each group. The figure indicates no no- 
ticeable systematic redshift dependence, except a 
hint that the high redshift slope is slightly steeper 
than the values at lower redshift in lest-V for the 
blue group, though the error is large. 

The offset nl^^^ is plotted in Figure 8 as a func- 
tion of redshift for each sample. Since the SBs 

depend on redshift for a given mass, we fit Attest 
Figure 8 with redshift dependence, 

txli, = MO + 2.5log{l + zT. (5) 

The results are shown in Figure 8 and Table 2. 

The results, corrected for dust extinction, are 
also listed in Table 2, where the extinction is given 
from the best fit SED model by K09. It should be 
noted that the dust extinction in galaxies is a less 
reliable parameter in SED fitting (K09). 

5. Discussion 

We have obtained the linear correlation be- 
tween the SB and stellar mass of galaxies in the 
GOODS-N region at 0.3 < 2; < 3, which is an 
analogy of the Kormendy relation. The correlation 
has a nearly constant slope independent of redshift 
and the color of galaxies, especially in the rest- 2; 
frame. Although the low-mass limit of the obser- 
vations depends on redshift, the similar slopes of 
the linear regression in different redshift bins are 
suggestive of a universal relation between the SB 
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Table 1 

Linear fit of the surface brightness in rest-V and rest-z frames to stellar mass of 

galaxies. 



rest-V rest-^; 
Redshift Sample N a Mrcst ^ ^ Mrcst 

(mag arcscc"^) (mag arcscc"^) 



0.25 < z < 0.50 


all 


491 


-1.16 ± 0.03 




U -V <0 


68 


-1.58 ± 0.09 




U -V >0 


423 


-1.21 ± 0.03 




M, > 10^" Mq" 


92 




0.50 < z < 0.75 


all 


626 


-1.08 ± 0.03 




U -V <0 


83 


-1.58 ±0.09 




U -V >0 


543 


-1.12 ± 0.03 




M, > 10^°Mq=- 


144 




0.75 < 2 < 1.0 


all 


930 


-0.99 ±0.02 




U -V <0 


152 


-1.38 ±0.08 




U -V >0 


778 


-1.02 ± 0.03 




M, > IO^^Mq" 


302 




1.0 < z < 1.25 


all 


698 


-1.07 ± 0.03 




U -V <0 


184 


-1.47 ±0.10 




U -V >0 


514 


-1.13 ±0.03 




M, > 10^° Mq" 


201 




1.25 < 2 < 1.5 


all 


396 


-0.99 ± 0.04 




U -V <0 


151 


-1.34 ±0.11 




U -V >0 


245 


-1.03 ± 0.06 




M, > IO^'Mq" 


124 




1.5 < 2 < 2.0 


all 


430 


-0.92 ± 0.05 




U -V <0 


174 


-1.31 ± 0.11 




U -V >0 


256 


-0.99 ± 0.06 




M, > 10^° Mq'^ 


177 




2.0 < z < 2.5 


all 


565 


-0.97 ±0.04 




U -V <0 


294 


-1.42 ± 0.07 




U -V >0 


271 


-1.09 ±0.06 




M, > 10^°Mq^ 


347 




2.5 < 2 < 3.0 


all 


184 


-1.04 ±0.08 




U -V <0 


127 


-1.21 ± 0.12 




U -V >0 


57 


-1.43 ±0.15 




M, > 10^" Mq^ 


104 





22 


.84 ± 0, 


,04 


0, 


58 


-1.32 ± 0.02 


22 


,10 ± 0, 


,03 


0.48 


21 


.68 ± 0, 


,16 


0, 


,46 


-1.40 ± 0.08 


21 


,60 ± 0, 


,16 


0.45 


22 


.86 ± 0, 


,04 


0, 


55 


-1.38 ± 0.03 


22 


,11 ± 0, 


,03 


0.45 


22 


.23 ± 0, 


,04 


0, 


,41 




21 


,28 ± 0, 


,05 


0.44 


22, 


.74 ± 0, 


.03 


0, 


,52 


-1.35 ± 0.02 


22, 


,05 ± 0, 


,02 


0.43 


21 


.80 ± 0, 


.12 


0, 


.39 


-1.55 ±0.09 


21, 


.57 ±0, 


.12 


0.38 


22, 


.77 ± 0, 


.03 


0, 


.51 


-1.39 ±0.02 


22, 


.07 ± 0, 


.02 


0.42 


22 


.12 ± 0, 


.05 


0, 


.61 




21, 


.22 ± 0, 


.05 


0.58 


22, 


.53 ± 0, 


.02 


0, 


.51 


-1.32 ±0.02 


21, 


.90 ± 0, 


.01 


0.42 


21 


.98 ± 0, 


,09 


0, 


,43 


-1.33 ± 0.08 


21, 


,78 ± 0, 


,08 


0.41 


22 


.56 ± 0, 


,02 


0, 


,52 


-1.36 ± 0.02 


21, 


,91 ± 0, 


,02 


0.42 


21 


.96 ± 0, 


,03 


0, 


,61 


21.08 ± 0.04 




0.62 






22 


.30 ± 0, 


,02 


0, 


,49 


-1.44 ± 0.03 


21, 


,72 ± 0, 


,02 


0.44 


21 


.84 ± 0, 


,07 


0, 


,45 


-1.58 ± 0.10 


21, 


,52 ± 0, 


,08 


0.47 


22, 


.37 ± 0, 


,02 


0, 


,48 


-1.48 ±0.03 


21, 


.76 ± 0, 


.02 


0.42 


21 


.67 ±0, 


.04 


0, 


,63 




20, 


.83 ± 0, 


.04 


0.63 


22, 


.08 ± 0, 


.03 


0, 


.51 


-1.35 ±0.04 


21, 


.57 ±0, 


.02 


0.44 


21 


.78 ± 0, 


.07 


0, 


,42 


-1.18 ±0.11 


21, 


.65 ± 0, 


.07 


0.44 


22 


.14 ± 0, 


,03 


0, 


,54 


-1.43 ± 0.05 


21 


,60 ± 0, 


,03 


0.44 


21 


.65 ± 0, 


,06 


0, 


,62 




20, 


,90 ± 0, 


,06 


0.64 


21 


.82 ± 0, 


,02 


0, 


51 


-1.36 ± 0.04 


21 


,33 ± 0, 


,02 


0.47 


21 


.56 ± 0, 


,05 


0, 


,45 


-1.39 ± 0.12 


21 


,30 ± 0, 


,06 


0.51 


21 


.91 ± 0, 


,04 


0, 


,52 


-1.38 ± 0.06 


21 


,35 ± 0, 


,03 


0.45 


21 


.40 ± 0, 


,05 


0, 


,62 




20, 


,66 ± 0, 


,05 


0.64 


21 


.38 ± 0, 


,02 


0, 


,47 


-1.41 ± 0.04 


21, 


,07 ±0, 


,02 


0.44 


21 


.21 ± 0, 


.03 


0, 


.42 


-1.51 ±0.08 


21, 


.03 ± 0, 


.03 


0.45 


21 


.58 ± 0, 


.04 


0, 


.46 


-1.44 ±0.05 


21, 


.11 ± 0, 


.04 


0.44 


20 


.91 ± 0, 


.03 


0, 


.61 




20, 


.39 ± 0, 


.04 


0.71 


20, 


.98 ± 0, 


.04 


0, 


.48 


-1.40 ±0.09 


20, 


.70 ± 0, 


.04 


0.51 


20 


.88 ± 0, 


,04 


0, 


,44 


-1.37 ±0.15 


20, 


,68 ± 0, 


,05 


0.53 


21 


.36 ± 0, 


,10 


0, 


,48 


-1.59 ±0.15 


20, 


,84 ± 0, 


,09 


0.46 


20, 


.57 ±0, 


.06 


0, 


,64 




20, 


.14 ± 0, 


.07 


0.72 



"■Average for galaxies with M» > 10^° Mq 
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Table 2 

Redshift dependence of the surface brightness. 



Sample n /.lq 

(mag arcsec""'^) 



rest- V 






all 


-1.74 ± 0.13 


23.63 ± 0.12 


U 




V > 


-1.46 ± 0.07 


23.49 ± 0.07 


U 




V <0 


-0.78 ± 0.25 


22.28 ± 0.24 


M, 


> 




-1.53 ± 0.14 


22.92 ± 0.13 


rcst-V (extinction corrected) 






all 


-1.97 ± 0.10 


23.13 ± 0.09 


U 




V > 


-1.78 ± 0.10 


23.04 ± 0.09 


u 




V <0 


-1.48 ± 0.26 


22.47 ± 0.25 


M, 


> 




-1.85 ± 0.10 


22.24 ± 0.10 


rest-2 






all 


-1.30 ± 0.11 


22.71 ± 0.10 


U 




V > 


-1.22 ± 0.09 


22.68 ± 0.08 


U 




V <0 


-0.83 ± 0.22 


22.13 ± 0.21 


M, 


> 


10^° Mq 


-1.05 ± 0.09 


21.75 ± 0.09 






rest-z 


(extinction corrected) 






all 


-1.42 ± 0.10 


22.44 ± 0.09 


U 




V > 


-1.38 ± 0.09 


22.43 ± 0.09 


U 




V <0 


-1.21 ± 0.24 


22.24 ± 0.23 


M, 


> 


10^°Mq 


-1.22 ± 0.09 


21.39 ± 0.09 



-0.9 


- (a) 


rest— F 








^ -1.2 












-1.5 














-0.9 


h 

- (b) 


rest— z 


— I 


o all 

□ U-V>0 

A U-V<0 


a -1.2 












-1.5 










i 



1 2 



Z 

Fig. 7. — Slope a of the linear regression in Equa- 
tion 4 as a function of redshift for each group in 
(a) rest-y and (b) rest-z frames. 




12 3 



z 

Fig. 8. — Surface brightness (^Jost) M, = 
as a function of redshift. Diamonds are 
the average SB for Af» > 10^° A/0. Extinction is 
not corrected. The line is the fitting result for each 
group (Table 2). 
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and stellar mass of galaxies in a wide mass range. 
The stellar system of galaxies keeps the tight cor- 
relation over the cosmic time. 

Considering that the light in the rest-z band 
is dominated by old stars and therefore traces the 
stellar mass, we note that the result should demon- 
strate a universal relation between surface stellar- 
mass density and stellar mass of galaxies. This 
finding indicates that stellar mass in galaxies was 
build up over cosmic time in a similar manner, 
irrelevant to galaxy mass, as opposed to the sce- 
nario that massive galaxies mainly accumulated 
their stellar mass by major merging. Dissipation- 
less merging increases the size of galaxies, which 
leads to lower SB (Navarro 1990; Boylan-Kolchin, 
Ma, & Quataert 2006; Mcintosh et al. 2005; Dam- 
janov et al. 2009; Nipoti et al. 2009). If dry merg- 
ing is a more important process for more massive 
galaxies, as conjectured from ACDM model, then 
the slope should become shallower at lower red- 
shifts for massive galaxies. 

Figure 8 shows the redshift dependence of /xj^g^ . 
As a whole, SB is brighter at higher redshifts for 
a given mass. The brightening obeys the relation 
with n = —2.0 ~ —0.8 in rest-V^ and n = —1.4 ~ 
—0.8 in rest-z, depending on the subgroups and 
extinction correction used. In Figure 9, we com- 
pare the observed SB of galaxies with the lumi- 
nosity evolution, for a given stellar mass with no 
size change, for two extreme cases, i.e. a constant 
star formation (SF) model and an instantaneous- 
burst galaxy model formed at Zf = 9 (Bruzual 
& Chariot 2003). It is noted that we discuss the 
properties of galaxies with same stellar mass at 
any redshifts, not tracing the evolution of mass 
or luminosity for a certain galaxy. The SB evolu- 
tion of red ([/ - F > 0) and blue {U ~ V < 0) 
samples is well reproduced by the luminosity evo- 
lution models of burst and constant SF models, 
respectively, though SB is a little bluer than the 
models at higher redshift. However, the models 
depend on various parameters (e.g., Zf , metallicity, 
star- formation rate, extinction, etc). For example, 
lower metallicity at higher redshift makes the color 
bluer, which gives an evolution more consistent 
with the observation. In any case, our result sup- 
ports the scenario that galaxies accumulate stellar 
mass mainly by star formation. 

The different selection criteria and analysis 
would not facilitate the direct comparison of the 
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Fig. 9. — Comparison of the observed surface 
brightness for blue, red, and massive samples with 
the luminosity evolution, for a given stellar mass, 
with a constant star-formation model (dash-dot 
line) and an instantaneous-burst galaxy model 
(solid and dash lines) formed at Zr = 9 (Bruzual 
& Chariot 2003). The models are arbitrarily offset 
at the rest-y band to fit the observations. 
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present result with those of previous studies, which 
discussed luminosity evolution mainly in half-light 
radius based on images with a high spatial resolu- 
tion. The studies at high redshifts could be biased 
toward compact and massive galaxies. In contrast, 
we focus on more general galaxies and studied the 
evolution of their SB based on the Kron radius 
for a given mass, making use of our deep imag- 
ing data, which are favorable for diffuse galaxies. 
Despite of the difference, the present results are 
consistent with previous studies. 

For example, our result for the blue population 
supports the finding of Bardcn ct al. (2005) that 
the average surface brightness of disk galaxies in- 
creases with redshift by about 1 mag to z ~ 1 
for galaxies. Ferguson et al. (2004) and Akiyama 
et al. (2008) showed the size evolution of LBG 
galaxies at high redshifts, which indicates much 
stronger SB brightening than our results. Al- 
though the blue group includes most LBGs in 
our catalog (Ichikawa et al. 2007), our differ- 
ent selection criteria sample various populations of 
more general star-forming galaxies with low stel- 
lar mass, which would have more moderate evo- 
lution (Trujillo et al. 2004; Trujillo et al. 2006). 
The results for red massive groups are also in good 
agreement with those in previous studies. Early- 
type galaxies were brighter by 1 ~ 2 mag at z 1 
(Mcintosh et al. 2005) and z ~ 1.5 (Damjanov 
et al. 2009), which would correspond to our red 
group. 

Trujillo et al. (2007) for galaxies with M* > 
1O"M0 and Toft et al. (2009) for those with 
> 5 X 10i°Mo showed that massive galax- 
ies were much smaller in the past than their lo- 
cal massive counterparts and the evolution was 
particularly strong for the highly concentrated ob- 
jects, which are supposed to be high SB galaxies. 
Since our sample includes few such massive galax- 
ies, we calculated the average SB and its evolu- 
tion in galaxies with Af, > lO^^M© for compar- 
ison. The result is in good agreement with their 
results for resolved galaxies. Their finding that the 
evolution is stronger for spheroid-like or quiescent 
galaxies would be supported by our result for the 
galaxies with the dust extinction corrected (Table 
2). 

Finally we emphasize that the present study has 
not only confirmed the previous studies on the 
SB evolution of galaxies, but also demonstrated 



that the brightening dates back to z ~ 3. We 
should also mention about the point-like or un- 
resolved galaxies (r < 0".4), which are excluded 
from the analysis for the present stellar-mass Ko- 
rmendy relation. For example, early-type galax- 
ies (Damjanov et al. 2009) and compact galaxies 
(van Dokkum et al. 2008) tended to be a factor 
2-5 smaller in half light radius than local counter- 
parts of similar mass. The SBs of galaxies at the 
faint-end limit of each redshift bin in Figure 6 are 
seriously affected by background noise, while more 
massive galaxies would not. The real SB of the 
massive galaxies, which are well above the SB limit 
but unresolved in the present observation, could 
be much brighter. The samples of e.g., Akiyama 
et al. (2008), Trujillo et al. (2007), and Muzzin 
et al. (2009) would be such galaxies as those we 
did not resolve with our poorer image size. If they 
are observed with higher spatial resolution, they 
could be plotted further below the regression lines 
in Figures 4 and 5. Therefore, the stellar-mass 
Kormendy relation presented in this paper would 
rather show a lower boundary of the correlation of 
stellar mass and SB of galaxies; lower SB galaxies 
are prohibited over the boundary. 

This work has been supported in part by a 
Grant-in- Aid for Scientific Research (21244012) of 
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